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ABSTRACT 

We present an analysis of the faint M star population seen as foreground contami- 
nants in deep extragalactic surveys. We use space-based data to separate such stars 
from high redshift galaxies in a publically- available dataset, and consider the pho- 
tometric properties of the resulting sample in the optical and infrared. The inferred 
distances place these stars well beyond the scale height of the thick disk. We find 
strong similarities between this faint sample (reaching i'^^g = 25) and the brighter 
disk M dwarf population studied by other authors. The optical-infrared properties of 
the bulk of our sources spanning 6000A-4.5/im are consistent with those 5-10 mag- 
nitudes brighter. We also present deep spectroscopy of faint M dwarf stars reaching 
continuum limits of « 26, and measure absorption line strengths in the CaH2 and 
Ti05 bands. Both photometrically and spectroscopically, our sources are consistent 
with metallicities as low as a tenth solar: metal-rich compared with halo stars at sim- 
ilar heliocentric distances. We comment on the possible MACHO identification of M 
stars at faint magnitudes. 
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1 INTRODUCTION 

The advent in recent years of extremely deep imaging sur- 
veys, designed to probe distant galaxies through multi- 
waveleng th imaging extending deep into the infrared (e.g. 
GOODS. iGiavalisco et aLll2004 ). has made possible studies 
of extremely faint, and often extremely red, sources. In or- 
der to target high redshift galaxies, colour cuts have been 
developed to exclude the majority of Galactic stars - in- 
cluding the low metallicity dwarfs expected in the halo of 
our own galaxy. Nonetheless, it has become clear that an 
unexpected number of extremely red stars remain at faint 
magnitudes. The spectroscopic component of high-redshift 
galaxy searches have confirmed the nature of such stars, in 
several ca ses obtaining deep spectroscopy as well as ph otom- 
etry (e.g. IStanwav et all l2004al : IVanzella etHI bOOel ). but 
the overall population has not been systematically charac- 
terised. 

Low mass stars of class M or later comprise more than 
80 per cent of the stellar population of our galaxy. Faint 
stars of late-M, L and T classes probe the tail of the stel- 
lar mass function, bridging the divide between the main se- 
quence and low mass brown dwarfs. Such stars have unusual 
colours, arising from the presence of deep molecular absorp- 
tion bands in their cool atmospheres. 
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The faintness of these stars has made their investiga- 
tion challenging, particularly with increasing distance from 
our sun. Efforts to determine the properties of such stars in 
the halo and thick disk have concentrated on nearby stars 
with high proper motions indicative of a halo origin (e.g. 
iGizis fc Reidl l 1999*1. While deep surveys of the M star popu- 
lation in order to identify distant examples have been under- 
taken with WFPC2 on the Hubble Sp ace Telescope (HST), 
reaching /as=23.2 (|Zheng et aLlbOOll ). these have been lim- 
ited to wavebands shortwards of 9000A where the coolest 
stars have very little flux. 

Meanwhile, the more local population - late M stars at 
disk metallicities within 1.5 kpc of the sun - have been the 
subject of detailed study in rccentyears. The Sloan Digital 
Sky Survey (SDSS, York et al. 20doi) is relatively shallow in 
the optical bands, reaching an / band limit of i' ~ 20. How- 
ever, the addition of a redder z' band in the optical to sim- 
ilar depths has enabled the photometric selection of cooler, 
class L and T stars as sources with large flux decrements be- 
tween adjacent passbands, caused by sharply-deflned, strong 
absorption and transmission features. The z'-band imaging 
also allows the characterisation of late M stars by their i' — z' 
colour. Such surveys have enabled templates based on thou- 
sands of local (mainly disk) dwarf stars t o be constructed 
(jHawlev et al.ll2002l : ISochanski et al.ll2007l 'l. 

At the same time, advancements in infrared detector 
technology have made it increasingly easy to study these 
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faint, cool stars at wavelengths longwards of the optical. 
Their near and mid-infrared bands are dominated by molec- 
ular absorptio n, leading to extreme colours in infrared bands 
l|Patten et al.L 2006). Again, such studies have been limited 
to local stars, in order to attain the maximum possible signal 
to noise. 

However, given the lower typical metallicity of stars in 
the halo when compared to the galactic disc, the character- 
istics of local disk-dominated M star samples may well be 
unrepresentative of fainter sources at larger distances and 
the colour selection methods applied locally might be ex- 
pected to detect nothing at halo distances. 

In this paper we present an analysis of faint low mass 
stars, selected for their extreme colours between two adja- 
cent bands in high spatial resolution, deep optical imaging. 
Such a colour selection is expected to select either relatively- 
nearby late M stars at sub-solar metallicities or near-solar 
metallicity early M stars at the large distances more nor- 
mally associated with the Galactic halo, if such a popula- 
tion exists. These two alternatives cannot be distinguished 
on the basis of a single red colour, but can be explored using 
spectroscopy and photometry probing well into the infrared. 

In section [5] we present the data used in this analysis 
and we discuss the selection of stellar candidates in section 
[3l In sections |4] and [5] we discuss the optical and infrared 
properties of our sample, and in section |S] we examine the 
photometry of cooler stars in the same fields. In section [7] 
we discuss existing spectroscopy of such faint stars, includ- 
ing a spectroscopic sample presented here for the first time, 
and spectroscopic abundance indicators that may hint at a 
typical metallicity for this population. Finally, in section [4] 
we discuss the interpretation of the overall properties of the 
faint, cool star population including their metallicity and 
distance distribution. 

All magnitudes in thi s paper (optical and infrared) are 
quoted in the AB system l|Oke fc Gunnlligsal ). 



2 DATASETS USED IN THIS ANALYSIS 

2.1 Filters and Photometric Systems 

We choose to use the AB system (|Oke fc GunnlligSSl ) when 
expressing magnitudes in this analysis. The reason for this is 
twofold. Firstly the AB system is widely used by extragalac- 
tic astronomy and the surveys that generate much of the 
deepest optical imaging, including the data discussed below. 
The AB magnitude scheme is also used in more local large 
spectroscopic samples such as the Sloan Digital Sky Survey 
(SDSS). Secondly, AB magnitudes are also constructed on a 
physical basis. A source flat in flux as a function of frequency 
(fi^ oc A~^) has zero colour in all bands. 

Colours are calculated throughout this paper by con- 
volving available observational and theoretical spectra with 
the filter transmissions concerned (including the instru- 
ment and detector response) to accurately account for 
filter-specific effective wavelengths and full-width half- 
maxima. Our main reference spe ctra for faint st a rs the 
observationally-derived spectra of iBochanski et all (|2007l ) 
which were constructed for each subtype of class M and later 
stars from a total sample of over 700 observed spectra taken 
by the SDSS. As figure [2] illustrates, these template spectra 
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Figure 1. The transmission profiles of the {v, dotted), 

F7To W (i' , dashed ) and F850LP {z' , dot-dash) wavebands, and 
of the lBesselll l ll990l) V (dark grey) and / (light grey) filters. Offset 
and plotted below is the spectrum of an M4 star. 

trace the same stellar locus as those of iPickled l| 19981 ). al- 
though slight deviations are seen between the two template 
sets, particularly at late spectral types. 

Use of the AB (rather than Vega) magnitude system 
introduces a blueward shift of 0.43 magnitudes in the "Stan- 
dard" Bessell V — I colour, which accounts for the spectral 
slope of a Lyr at long wavelengths. In addition the F&Q&W 
(v) and F775W (i') filter set has redder central wavelengths 
than their Bessell filter counterparts, and the FQQ&W filter 
is broader than the "Standard" V as figure [T] illustrates. 

2.2 Surveys 

The analysis presented in this paper utilises two deep, multi- 
wavelength imaging surveys in order to combine their differ- 
ent strengths. 

The photometric selection of cool dwarf stars paral- 
lels that of high redshift galaxies in many respects. In both 
cases the accurate characterisation of the sources, and their 
reliable separation from other sources with similar optical 
colours, requires multiwavelength imaging spanning a long 
baseline from 5000A to 3.6/im, and high angular resolution 
impossible to obtain from the ground without the use of 
adaptive optics. 

The Great Observatories Origins Deep SurvejQ 
(GOODS) satisfies the requirements for such a selection. In 
addition to an abundance of auxiliary data, the survey com- 
prises two primary data sets, optical and infrared, each of 
which surveys a field in the northern hemisphere (GOODS- 
N) and a second in the south (GOODS-S). AU GOODS data 
is publically released to the astronomical community in the 
form of fully reduced images. 

The optica l (0. 3-l^tm) component of GOODS 
iGiavalisco et al.1 |2004| ) consists of deep optical imaging 

^ http:/ /www. stsci.edu/science/goods/ 
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in the F435W^ (b), {v), F775W {i') and F850LP 

{z') bands, each reaching depths > 28th magnitude as 
shown in table [T] obtained using the wide field channel of 
the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). The imaging has been drizzled to a 
pixel scale of 0."03 per pixel, with a point source FWHM 
of 0."05 in the i' band. In addition, the GOODS team 
have made publically available catalogues of the source 
photometry and properties in ea ch field, generated usin g 
the SExtractor software package (|Bertin fc Arnout3ll996h . 
and with detection parameters tuned to the survey depth 
and resolution. In this analysis we utilise the Kron-radius 
based automatic magnitudes reported by SExtractor 
(i.e. MAG_AUTO). These magnitudes follow the curve 
of growth of the source profile, and are equivalent to a 
corrected aperture magnitude for an unresolved source. We 
use version rl.l of the GOODS catalogue^. 

The second primary dataset generated by the Great 
Observatories Origins Deep Survey comprises deep Spitzer 
Space Telescope in frared (3-10/im) imaging of the fiel ds tar- 
geted by Hubble (|Dickinson fc GOODS TeamI |2005| ). The 
IRAC instrument was used to survey the GOODS fields 
at 3.6, 4.5, 5.8 and 8.0 /xm, with a pixel scale of 1.2"/pix 
to a depth >26th magnitude in the first two bands and 
>24th magnitude at longer wavelengths as shown in table 
As discussed in section (5] confusion is a significant is- 
sue for sources at the faint magnitudes of our targets, and 
in imaging of this depth. As a result, the 12 arcsecond 
(10 pixel) apertures recommended for IRAC photometry in 
spars e fields seldom esc ape confusion in the GOODS im- 
ages. IVerma et al.l (|2007l ) determined that smaller apertures 
with a diameter of 4.5 arcseconds were suitable for compact 
sources in the GOODS images, and could be reliably cor- 
rected for fiux in the wings of the point spread function. 
We apply th e same prescription for aperture magnitudes as 
IVerma et al.l . correcting to total magnitudes using the off- 
sets given in table [T] 

Both GOODS fields look out of the plane of the galactic 
disc and away from the galactic centre, with galactic coordi- 
nates of 1= 125.865, fe=-H54.808 for the GOODS-N field and 
^=223. 57, 6=-54.43 for the GOODS-S. Both fields were also 
selected to lie in regions of low Galactic extinction, selecting 
against structures in the galactic disk. 

The spectroscopic analysis of faint M stars presented in 
section[7]is also based in part on GOODS data. Several of the 
stars identified in section [3] have deep optical spectroscopy 
taken with F0RS2 at the Very Large Telescope (VL T) as 
part of the ESO GOODS survey (|Vanzella et al.ll2006D . 

In order to improve the statistical significance of our 
analysis, we supplement the ESO GOODS data with fur- 
ther M class star spectra observed as part of the BDF 
project. The BDF survey comprises deep multicolour imag- 
ing and spectroscopic follow-up of red sources in four near- 
contiguous fields, applying an 7? — / colour selection to iden- 
tify galaxies at 2 > 5, and reaching a limiting depth of 
Tab = 26.3. The fir s t of th ese fields, BDFl, was described in 
iLehnert fc Bremed l|2003D while the remaining three fields, 
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^ The GOODS programme also obtained imaging with the MIPS 
instrument that is not discussed here 



BDF2-4, were observed with an identical strategy and will 
be described in a forthcoming paper. These data are not used 
in the photometric analysis since with ground-based seeing 
it is impossible to separate mid-M stars from high redshift 
galaxies on the basis of optical colours alone (without B 
band imaging reaching some four magnitudes deeper than 
the / band limit). However, the high spectroscopic complete- 
ness of this survey for objects with extreme R — I colours 
(essentially the same selection function as the HST/ ACS 
V — I selection discussed in section [Sjl makes it extremely 
useful for the a posteriori analysis of stellar spectra. 

The FORS2 R (R_SPECIAL, ESO number 76) and I 
(IJ3ESS, ESO number 77) band filters used for imaging in 
the BDF fields were selected to be sharp sided and to have 
minimal overlap. Despite the different naming, the wave- 
length coverage of the R band at 6550A (FWHM=1650A) 
overlaps significantly with that of the v band at 5907A 
(FWHM=2343A). Hence a colour cut of (i? - I)ab > 1.5 
reproduces the HST/ ACS v — i' > 1.3 selection function for 
stars discussed in section|3] selecting class M3-M4 and later. 



3 A PHOTOMETRIC SAMPLE OF FAINT 
M-DWARFS 

Class M, L and T stars have an underlying red continuum 
due to their cool temperatures. Their spectrum in the yellow 
and red is dominated by jagged molecular absorption bands 
that further redden the v — i' and V — I colours. As a result, 
such cool stars can have extreme colours in adjacent filters 
and are regularly selected by 'dropout' surveys searching for 
high redshift sources (which show a similar drop between 
bands due to neutral hydrogen absorption in the intergalac- 
tic medium). We note that although 'dropout' is a historical 
term used by convention in the selection of high redshift 
galaxies, neither high redshift galaxies nor cool stars are ex- 
pected to have zero flux in the bluewards band. As a result, 
whether the sources show a finite colour or drop entirely be- 
low the detection limit in the bluewards band depends on 
the relative depth of the imaging. Hence the term 'break' (as 
in Lyman break galaxy) or 'drop' is more technically correct. 
Sources are required to satisfy the criterion of a significant 
spectral break, measured through their photometry, but not 
to remain undetected in any given band. Nonetheless, for 
historical reasons, we continue to use the term 'dropout' 
while cautioning readers that its literal meaning is inaccu- 
rate. 

In order to select mid-M class stars and later, we apply 
the same u-drop selection function commonly used to select 
high redshift galaxies. As figure [5] illustrates, a colour cut 
of {v — i')AB > 1.7 effectively selects stars of class M3 
and later at solar metallicity (see section |4] for discussion 
for extreme low metallicity populations which will fall out 
of our colour selection). We note that we are likely to be 
insensitive to later spectral types: despite their red colours 
in V — i' , these sources are generally too faint in i' to be 
detected to our survey limit in the i'-band. Beyond M6, 
the almost-linear increase in v — i' colour begins to plateau, 

which has been used to identify the Lyman-a spectral break in 
galaxies at 4.8 < z < 5.8 
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Band b v i' z' 3.6 /^m 4.5 /xm 5.8 /im 8.0 

Depth 28.87 29.46 28.85 28.55 26.96 26.38 24.50 24.37 

mtot-map - - - - -0.180 -0.180 -0.326 -0.418 

Table 1. The 3(t depth of the GOODS imaging used in this analysis, and aperture corrections appHed in the Spitzer wavebands. The 
Hmits for optical bands are given in a 1 arcsecond aperture and those for the Spitzer bands in a 4.5 arcsecond aperture, corrected to 
total magnitude. 
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Figure 2. The v — i' and i' — z' colours of cool (class M and L) 
stars, calculated from empirical template spectra and convolved 
with the HST/ AC S filter set. So lid points indicate the cool star 
sample of iBochanski et al. I l|2003), diamonds indicate the M star 
spectra of lPickle^Hl99^ while crosses are K stars from the same 
source. M class stars increase linearly in i' — z' colour with in- 
creasing subclass, although we note that the two sets of empirical 
spectra differ significantly at certain subclasses. 



reducing the ability of this colour to distinguish between 
M class subtypes. Adding a second, redwards colour such as 
i' ~z' improves the discrimination between cool stars. Hence 
in order to examine the coolest stars in our survey fields we 
also consider a z'-band selected sample of stars satisfying 

{i'-z')AB > 1.30. 

The resulting sample is expected to contain several dis- 
tinct populations: the targeted stars, high redshift {z > 
5) galaxies and quasars which show a strong break at 
Arest =1216A, and old elliptical galaxies at 1 < z < 2 which 
break at Xrest =4000A. The majority of contaminant galax- 
ies can be identified by their extended half hght radii in 
the HST/ ACS data. All confirmed z > 5 galaxies selected 
by such a colour cut criterion are reso lved in space-based 
imaging, albeit barely in several cases jBunker et aLll2003l: 
Bremer et al. I l2004l : IStanwav et al. I l2004al lbl: iBouwens et al.l 
20041 ). Similarly, galaxies at 1 < 2 < 2 are typically well 
resolved by HST. As figure |3] illustrates, unresolved sources 
separate cleanly from extended sources in the GOODS imag- 



GOODS-S 
GOODS-N 



150 
150 



53 
71 



Table 2. The number of unresolved sources in the deep 
HST/ ACS imaging satisfying our colour selection criteria as de- 
scribed in section |3] 



ing for sources brighter than = 25 (or for — 26 in 
the case of z'-selected sources) . 

Clearly high redshift quasars cannot be distinguished 
from stars by their half hght radii, since both populations are 
expected to be unresolved. The predicted number of quasars 
at the faint magnitud es probed here is expected to be small. 
iDouglas et al.l (|2007f ) identified only one active galactic nu- 
cleus (AGN) in a highly spectroscopically complete survey 
of dropout objects examining 450 arcmin^ to a depth of 
26th magnitude in That source was slightly resolved 
in HST/ ACS data, with comparable contributions to the 
rest-frame ultraviolet fiux from the AGN and a starburst in 
the host galaxy. Hence we might expect to detect one AGN 
in the 300 arcmin^ GOODS survey area, either unresolved 
or marginally resolved. Fortunately, such a source should 
be di stinguishable based on its infrared colours (Ster n et al.l 
I2OO6I ) as discussed in section [5] 

In order to select low mass stars in the available ACS 
imaging with minimal contamination from extragalactic 
sources, we apply an additional selection criterion, requir- 
ing that: FWHM,/ 4.0 pixels (0".12), and we truncate 
our u-drop selection at i'j^g = 25 and our i'-drop selection 
at Zj^j^ — 26. 

Each source was visually inspected to ensure that the 
photometry was not contaminated by neighbouring objects. 
In this deep imaging, the area of sky affected by sources 
above the detection limit is 2%. 

The resultant number counts of dropout-selected stars 
is shown in table (2] In this analysis we focus on the M stars 
in the sample since these are both numerous, allowing a sta- 
tistical analysis, and relatively bright, allowing spectroscopic 
follow-up. 

While the i'-drop selected sample is too small to derive 
reliable statistics, the magnitude distribution of the u-drop 
sample (i.e. M stars) is illustrated in figure [l] There is no 
clear trend in the number density of cool stars with apparent 
magnitude. The northern GOODS field is 34% more abun- 
dant in low mass stars than the southern field, despite their 
similar areas and orientation with respect to the galactic 
disk. 



^ which identifies galaxies at 5.6 < 2 < 6.5 or class L and T stars 



Similar to the surface densities for faint AGN found by 



iMahabal et al.l ||2005| ') and lShankar fc Mathu"il l[20o3) 
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Figure 3. The distribution of Gaussian FWHM determined for 
sources in tfie GOODS fields. The stellar locus is clearly separated 
from the distribution of extended galaxies for i'j^g < 25 with very 
few ambiguous sources. We apply a criterion of FWHM < 4 pixels 
to select stars. Note that sources saturate in the GOODS imaging 
at 19th magnitude 
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Figure 4. The apparent magnitude distribution of v-drop stars in 
the GOODS field. The decline in number counts at bright magni- 
tudes arises both due to the small volumes probed at bright mag- 
nitudes and due to the potential saturation of brighter sources 
in the deep GOODS imaging. There is no obvious decline in the 
number of sources seen at faint magnitudes. Each field surveys an 
area of 150 arcmin-^. 



4 THE OPTICAL PROPERTIES OF FAINT 
M-DWARFS 

The selection of stars by photometric methods relies on a 
combination of the red optical spectrum that arises from a 
cool blackbody with the abrupt absorption edges of metal 



molecules. In the absence of significant metal enrichment, 
the colours of stars at a given temperature are likely to 
be less extreme since the strength of molecular features is 
reduced. Hence a sample of faint stars based on extreme 
colours in two adjacent bands is expected to select either 
relatively-nearby late M stars at low metallicities as well as 
near-solar metallicity early M stars at the large distances 
more normally associated with the Galactic halo, if such a 

populati on exists. 

The lAUard fc Hauschildtl lll995h 'NexGen' models and 
the Kurucz ATLAS9 models l|CasteUi fc Kurucj |2004| : 
iKuruc j \W9^ remain the most recent stellar atmosphere 
models to attempt the temperature and metallicity regime 
of interest here (although the ATLAS9 models in fact do 
not survey T<3500K late M stars and do not account for 
the important Ca H bands in cool stel lar spectra). As the 
recent analysis bv lBertone et all |2004| ) shows, neither pro- 
vides a satisfactory fit to the spectra of cool stars. While the 
MARCS spectra produce a somewhat better fit, the publi- 
cally available library only extend as cool as 4000K. In prac- 
tise, no stellar models can reliably reproduce the colours of 
cool stars. 

On the other hand, there exist a number of published 
'extreme' cool stars known to be at significantly sub-solar 
meta llicity (e.g. lGizi^ll997l : [Burgasser. Cruz, fc Kirkpatrickl 
'2007*). Unfortunately the photometry and/or spectral cov- 
erage of these observations inevitably covers too short a 
wavelength range to reliably calculate v, i' and z' colours 
simultaneously from the empirical spectra. 

As a result, while we tentatively consider the results 
of available model templates as well as empirical spectra, 
we caution the reader that such results should not be con- 
sidered reliable. We note that future interpretations of this 
population would benefit from improved modelling and in- 
terpretation of these faint stars in the red spectral region 
(>8000A) accessible to red-sensitive spectrographs. 

We have determined the optical colours expected for 
cool stars, in the same filter set employed by the GOODS 
survey, over a range of me tallicities. We use the M-star spec- 
tral models of lAllard fc Hauschildt (1995|!| over the range 
[Fc/H]=0. to -4.0 (where such models are available) and at 
four different temperatures: T=2000K, 2500K, 3000K and 
3500K. 

Figure [5] illustrates the colours predicted by the 
'NextGen' models, and compares these wi th the colours of 
the SDSS dwarf star templates derived bv lBochanski et al.l 
(|2007) ) which are dominated by nearby disk stars and the 
photometric sample derived in section O It is clear that 
while the v — i' colour fails to distinguish between differ- 
ent metallicity populations, the addition of an i' — z colour 
may provide a crude measure of metallicity. We note that 
a slight offset (approx 0.1 magnitudes in i! — z') appears 
to be required to match the model spectra to the observed 
near-solar, old disk stellar sample at bright magnitudes. 

Given the uncertainties in current models of M star at- 
mospheres, we compare their predictions with observational 
data on Galactic globular clusters. Observational evidence 
for the i' — z' (or more generally I — Z) colours of cool stars 



^ In their 'NextGen' form as implemented by 
iLeieune. Cuisinier. fc Bused lll998ll 
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Figure 5. The v — i' and i' — z' colours of model dwarf star atmospheres (from lAUard &: Hauschildtlll995l . , triangles), convolved with 
the HST/ ACS filter set. Points at different temperature for the same metallicity are joined by lines and the temper ature in kelvin is 
labelle d. Figure a (left) overlays these tracks with the colours of moderate metallicity dwarf star templates derived bv lBochanski et al.l 
while figure b (right) overlays the colours of faint stars identified in section |3] Photometric errors in both colours are <0.05mag. 



as a function of metallicity seems to be poor. However the 
V — I colours of such stars (observed in the Johnson-Cousins 
or Bessell systems and in Vega magnitudes) have been re- 
ported from observations of several globular clusters. While 
few observations probe deep enough to identify the base of 
the main sequence (as opposed to the main sequence turn- 
off), some observations of low metallicity globular clusters 
have characterised the main sequence ridge to very faint 
magnitudes. 

NGC 6397 is one such cluster, with a metallic- 
ity of [Fe/H]=-2.0 and ly ing at a distance of 2.2 kpc 
l|Cool. Piotto. fc Kind [1999 ). The main sequence in this 
source has been observed with HST/ ACS down to a limit- 
ing magnitude of /f814W ~ 30 (Vega) and is seen extending 
redwards to Vpeoew —Ifsuw ~ 4, the theoretical hy drogen- 
burning limit at this metallicity (|Richer et al.ll2006l ). At our 
colour cut of ti — i'=1.7 (approximately V — /voga ~ 2.6), 
the main sequence in this cluster corresponds to an ob- 
served magnitude of approximately /_814 — 24.5 (Vega), 
or an absolute magnitude limit of 12.4 in the I band, while 
our detection limit of i'j^g = 25 would only select sources 
with V — /vega < 3.0. This cutoff in apparent magnitude 
implies that we are only sensitive to approximately class 
M4-M5 stars at the distance and metaUicity of NGC 6397 
l|lfawlev et al.|[20o3 ) and would not expect to select any stars 
with halo metallicities at greater distances. This emphasises 
that our colour cut allows only for a narrow range of low 
metallicity objects which are sufficiently red and sufficiently 
bright to meet our colour selection criteria. 

In figure [5] we illustrate the effects of decreasing metal- 
licity on the main sequence colours of globular clusters with 
published photometry. As discussed above, NGC 6397 is a 
distant globular cluster with the low metallicity expected of 
halo stars. By contrast NGC 6366 and NGC 6144 are two 
clusters with [Fe/H]=-0.58 and [Fe/H]=-1.7 respectively. 



and published main sequence star photometry to fa,int m ag- 
nitudes observed with HST/ ACS (ISaraiedini et ahllioOTl ) as 
part of an ACS survey of globular clusters. While these data 
are amongst the deepest available, they do not reach the red 
colours, and hence faint 't^-band magnitudes of the M stars 
of interest here. Nonetheless, it is interesting to note that 
the well-defined main sequence in these globular clusters 
marks a gradual shift to bluer V — I colours with decreasing 
metallicity, and all three are bluer than the old-disk stars 
measured by the SDSS. 

The spectra of several very low metallicity ex- 
treme subdwa rfs in the field have n ow been published. 
Burgasscr, Cr uz, fc Kirkpatrickl (|20o3) demonstrate the ef- 
fects of decreasing metallicity on the optical spectra of late 
M dwarfs. The weakening of metal oxide bands at 7200A, 
7900A and 8500A with decreasing metallicity boosts stellar 
flux in the i' and z' bands. In particular, the deep 7900A ab- 
sorption feature in the spectra of cool stars at disk metallici- 
ties suppresses the i'-band magnitudes of these sources. De- 
creasing metallicity removes this suppression, boosting the 
v — i' colour of these sources while decreasing their colour in 
i — z' . While none of the published extreme subdwarfs have 
publishe d z'-band photometry, synthetic colour s calculated 
from the iBurgasser. Cruz, fc Kirkpatrickl l|2007l ) spectra in- 
dicate that at these low metallicities, late M dwarfs are ap- 
proximately 0.2 mag bluer in i' — z' than typical old-disk M 
stars (Bessell, private communication). 

It is clear from the combination of theoretical and ob- 
served data that most stars earlier than M3-M4 will be omit- 
ted at low metallicities due to their blue colours, while later 
M class stars are only likely to be detected relatively nearby 
due to their faint magnitudes. Thus our selection preferen- 
tially identifies stars with disk metallicities, while probing 
to halo distances. This analysis of our colour selection func- 
tion is supported by the distribution of optical colours in 
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Figure 6. T he colour-absolute magnitude sequenc e of three glob- 
ular clusters l l Richer et al .12006 ':' Saraicdini ct al. 2007) and of the 
cool stars identified by the SDSS (Bochanski ct al. 2007). In each 
case information in the given references is used to correct sources 
to absolute magnitudes on the AB system and, where necessary, 
to correct for Galactic reddening. For ease of comparison with the 
globular cluster work, the SDSS stars are shifted into the HST 
ACS F606VK and F814V1/ bands. The data in NGC 6144 and NGC 
6366 become unreliable redwards of approximately V — / = 1.3 



our faint star sample as shown in figure[S}3. The majority of 
stars have colours consistent with metallicities in the range 
— 1 < [Fe/H] < (i.e. similar to that of the Galactic disk). 

Two stars in the GOODS-N field show anomalously low 
i! — z' colours (j! — z' < 0.5, v — i' = 2.1) as might be expected 
at [Fe/H]>-2, and both of these are relatively faint {i' > 23). 
Comparison with the colours o f the e xtreme subdwarfs of 
iBurgasser. Cruz, fc Kirkpatrickl (|2007| ) suggests that these 
might be identified as very low metallicity sources, although 
both are somewhat redder in v — i' than might be expected. 
Either spectroscopy of these two blue stars, or optical pho- 
tometry of known extreme subdwarfs will be required to 
make a more quantitative comparison. 



5 THE INFRARED PROPERTIES OF FAINT 
M-DWARFS 

Having identified the stars in HST data, we now consider 
the colours of these unresolved sources in the infrared with 
Spitzer GOODS data. The deep images are close to the con- 
fusion limit in the IRAC bands, and many of the w-drop 
sources are confused in the imaging. Of a total of 124 v- 
drop stars in the GOODS fields we classify 48 (39%, figure 
[7] case a) as isolated, a further 36 (29%, figure [7] case b) as 
confused but potentially deblendable and the remaining 40 
(32%, figure[7]case c) as hopelessly confused. In this analysis 
we use only isolated objects since the noise in deconvolved 
magnitudes at this level is often dominated by residuals from 
neighbouring objects rather than the target source. 

The infrared colours of the remaining sample are shown 
in figures [8] and |9] with the c olours measured for M , L and T 
dwarfs in the same bands bv IPatten" et al.l l|2006h shown for 
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Figure 8. The [4.5]-[5.8] versus [3.6]-[4.5] colours of isolated stars 
selected a s ri-drops in the IR AC GOODS fields. The cool star 
sample of iPatten et al.l ||2006|) which comprises M5-T9 stars at 
bright magnitudes is shown for reference. A typical uncertainty in 
the colours (dominated by uncertainty in the aperture correction 
and residual flux from nearby sources) is shown in the upper left. 



comparison. Working some three magnitudes above the 3 a 
background limit, the errors in the colours are dominated 
by uncertainty in the aperture correction and are set to an 
indicative 0.1 mag in each band. 

The colours of cool stars in the infrared are dominated 
by the relative strengths of molecular absorption bands in 
their atmospheres, in particular those of CH4 in the 3.6 
band, CO in the 4.5 fim band and H2O in the 5.8 fim band. 

The infrared colours of our faint v-drop selected sources 
ar e very similar to those of the brighter population studied 
bv IPatten et al.l . despite being fainter by some ten magni- 
tudes. The GOODS u-drop stars extend slightly blueward 
of the bright sample in the [4.5]-[5.8] band, as would be ex- 
pected given that our u-drop sample extends to earlier M 
subclasses than the bright sample (M3/M4 as opposed to 
M5). It is not clear from existing models or data whether 
a similar bluewards shift might be expected with varying 
metallicity. 

While the bulk of the population has colours consistent 
with brighter M stars, a subset of seven sources are anoma- 
lous with M star optical colours and T star colours in the 
Spitzer wavebands. These fall loosely into three categories. 

One source, labelled 'A' on figure has colours con- 
sistent with a late-M/early-L classification within approx- 
imately 3ct of the photometric errors in both optical and 
infrared bands. The colours are also consistent with those 
expected for a close M/L binary system with the M star 
flux dominating the optical bands and the fainter L star 
contributing to the flux longwards of 3 ^im. Binary M/L 
and M/T syste ms have been observed at bright magnitudes 
in th e past (e.g. iLeggett et al.|[20o3 : IBurgasser fc McElwainI 
I2OO6I ). This source is towards the faint end of our sample, 
with i'j^g — 24.12. A binary with separation of l".09 at a 
heliocentric distance of 35 pc (using the M/L binary pair 
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Figure 7. Examples of GOODS/ ACS selected sources in the GOODS/IRAC imaging. Some sources are isolated or have only much 
fainter neighbours (e.g. case a), while others are hopelessly confused in the IRAC bands (e.g. case c) with neighbours of comparable 
magnitude or brighter. Some sources are blended but not beyond hope of deconvolution (e.g case b). Boxes are 20" to a side. 
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Figure 9. The [3.6]-[8.0] versus [4.5]-[5.8] colours of isolated stars 
selected as n-drops in the IRAC GOODS fields. As in figure [8] 
Error bars are shown separately for stars with anomalous colours 
in these bands (i.e. [3.6]-[8.0] >-l). These fall into three categories; 
A, B and C as described in section [5] 



of iBurgasser fc McElwairJ[2006l . as a model) would be un- 
resolved at HST/ ACS resolution at the distance of Ikpc 
estimated for an M6 star at this magnitude. If it is indeed 
a binary system then each component will be intrinsically 
fainter still and may lie at greater distances than estimated. 
Radial velocity measurements of the system may be neces- 
sary to determine whether the star is an anomalous M dwarf 
or a binary pair. 

A further three sources, labelled 'B' on figure |9] have 
z'-5.8 colours consistent with those of M stars but are 
anomalously bright in the 8 /im band. In one case, the pho- 
tometric errors are large and the faint star is within 2 a of 
the colours expected for an M star. However, the other two 
stars are well detected in all bands, and show no evidence 
of unreliable photometry. Both the 3.6 and 8 fim bands are 
dominated by CH4 absorption with the fundamental band 



in the bluewards filter. There is no clear interpretation of an 
8/im excess. 

Finally the remaining three sources, labelled 'C on fig- 
ure (9] show no distinctive breaks or obvious features in their 
spectral energy distributions but their flux peaks further to 
the red than expected for M stars. Again these sources have 
reliable photometry with no clear reason to believe photo- 
metric errors are larger than shown. As such they have M 
star colours in the optical and T stars in the infrared. Inter- 
estingly their i' — [3.6] colours are well within the distribution 
of the rest of the sample, showing no evidence for M/T mul- 
tiplicity (which would lead to red colours as the T dwarf 
contribution becomes dominant at 3.6 fim). 

We note that all of these anomalous sources defy 
straightforward explanation based on photometry alone and 
may benefit from spectroscopic analysis. 

One unresolved object in the GOODS-N has extremely 
unusual colours in the IRAC wavebands, lying outside of 
the region plotted in figures |8] and |9] The colours of this 
source {{z'-[3.6])ab=0.71, ([3.6]-[4.5])Ai3=0.09) maybe con- 
sistent with identification as a high redshift quasar or com- 
pact galaxy and inconsistent with any cool star. 



6 THE PROPERTIES OF COOLER STARS IN 
THIS SAMPLE 

To a limit of z'^^g = 26, there are a total of 12 unresolved 
sources in the GOODS fields that are identified as i'-drops 
but do not form part of the u-drop sample due to their faint 
magnitudes in i' . 

As figure [2] indicated, an i'-drop selection is sensitive 
to stars of spectral classes later than M8 (assuming near- 
solar metallicities) . The linear increase of i' — z colour with 
spectral type becomes less clear cut at the same subclass, 
due to the difficulty in assigning spectral types based on 
spectroscopic indices, making classification of these sources 
impossible based on optical photometry alone. 

Seven of these sources are unconfused in the 3.6 /xm 
band, although given their faint magnitudes, one is un- 
detected at 4.5 /xm a further four undetected longwards 
of 4.5^111. Figure [TU] shows the [4.5]-[5.8] versus [3.6]-[4.5] 
colours or 3 (J limits on the colours of the six stars for which 
these colours can be determined. All six are consistent with 
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Figure 10. The [4.5]-[5.8] versus [3.6]-[4.5] colours (open circles) 
or 3 (T limits on the colours of isolated stars selected as i'-band 
dropouts. Given their i' — z' colours, these sources are expected 
to be of class M8 or later. Symbols as in figure [8] 



the infrared colours of local M and L dwarf stars studied by 
iPatten et all (|2006l ). 

We note that an i'-drop selection used to identify cool 
stars here is likely to be significantly incomplete for low 
metallicity dwarfs if the available stellar models provide a 
reasonable approximation to their colours, as figure [S] illus- 
trates. A number of low metallicity extreme subdwarfs have 
now been found. Unfortunately none of these have published 
z-band photometry. 

.Burgasser. Cruz, fc Kirkpatrickl (|2007l ) present the op- 
tical spectra and I — J colours of a sample of late M and L 
extreme subdwarfs. The spectroscopically confirmed subd- 
warfs show I — J colours between 1.4 and 1.8 for M dwarfs 
increasing to about 2.4 for L4 dwarfs. Given the spectra 
of M class stars are essentially flat in fA in the range 8000- 
lOOOOA, it seems unlikely that most late M subdwarfs would 
have i' — z' colours as extreme as 1.3. Indeed, synthetic 
colours calculated from these spectra indicate that at these 
low metallicities, late M dwarfs are approximately 0.2 mag 
bluer in i' — z than typical old-disk M stars (Bessell, pri- 
vate commu nication). Nonetheless, the early L class extreme 
subdwarfs of iBurgasser. Cruz, fc Kirkpatrici3 (|2007l ) lie well 
redwards of the colour selection criterion used here. 

This metallicity-dependent selection effect complicates 
the analysis of L and T star number counts, and larger sur- 
veys will be needed to accurately constrain the surface den- 
sity of these cool stars at faint magnitudes.. Further progress 
in understanding the incidence and properties of low metal- 
hcity class M, L and T dwarfs may be possible in future given 
improved modelling of the effects of varying metallicity on 
the spectrum longwards of 1 [xm and in the Spiteer/IRAC 
wavebands. 



T SPECTROSCOPY OF DISTANT M-DWARFS 

Some tens of faint low mass stars have now been observed 
in deep spectroscopy, as part of the follow up to dropout- 
selected high redshift surveys. 

Extremely deep 8 m class telescope spectra of M or L 
class star s with i! > 22 have bee n pu blished or publically re - 
leased bv lStanwav et all (Hooii) and lVanzella et al.l (|2006l ). 
Additional deep stellar spectra were ob tained as part of the 
ESQ Remote Galaxy Survey (ERGS, lOoug las et al' '20071. 
Douglas 2007, in prep) and bv lLehnert fc Bremer (2003.) as 
part of their BDF survey. Many such spectra remain unpub- 
lished, being considered incidental to the search for z > 5 
galaxies. 

Six of the v-drop sample discussed above have deep 
spectroscopy obtained in multi-object mode with F0RS2 on 
the VLT as part of the GOODS programme and publically 
released bv lVanzella et al] (|2006l ). These sources have mag- 
nitudes in the range i'j^g = 23 — 24.5 and were observed at 
a dispersion of 3.2A/pixel and a spectral resolution R=660, 
for a typical exposure time of 14400 seconds. 

As discussed in section[2] we supplement this small sam- 
ple with a further 18 faint M class star spectra observed 
as part of the BDF project. The BDF survey comprises 
deep multicolour imaging and spectroscopic follow-up of red 
sources in four near-contiguous flelds, applying an i?-drop 
colour selection to identify galaxies at z > 5 and including 
a number of stars (which can be difficult to separate from 
compact galaxies on the basis of optical colour alone). 

i?-drop sources were observed spectroscopically in 
multi-object mode. Two FORS2 slitmasks were observed in 
each field, with a total integration time of 13000 seconds 
per mask, divided into twenty exposures of 650 s, each off- 
set by a different distance along the slitlets. Each slitlet was 
1 arcsecond wide. The spectra were wavelength calibrated 
from arclamps and flux calibrated from standard star ob- 
servations taken at the same time. The flnal spectra span 
the spectral range 6000-lOOOOA at a dispersion of 3.2A, and 
with a spectral resolution R=6 60. This configu r ation is vir- 
tually identical to that used bv lVanzella et al] (|2006l ). with 
the two samples differing only in the selection filters. 

Examples of such faint star spectra are shown in figure 
111! As expected for sources a hundred times fainter than 
those observed by the SDSS, the spectra are noisy. However, 
despite the low signal to noise of such spectra, they show 
the similar molecular absorption features to their brighter 
counterparts. Each of the twenty-four stars in the final spec- 
troscopic sample discussed here has i'j^g > 23 and signal to 
noise on the continuum around the features of interest (i.e. 
around 7000A) of 10 or greater (in some cases higher than 
20). 

Given these deep spectra it is possible to identify a num- 
ber of molecular absorption features, both lines and absorp- 
tion bandheads. Here we consider two diagnostic features in 
the spectrum of M stars: the surface gravity dependent Na I 
doublet at 8183, 8195A and the metallicity dependent TiO 
to CaH ratio around 7000A. 

The Na I doublet is a strong feature in the spectrum of 
M stars, but its measurement is complicated by the crowded 
spectrum around 8200A wh ich makes the tr u e con tinuum 
level difficult to determine. iFaber fc Frenciil (|l980l ) intro- 
duced a systematic measurement of the equivalent width 
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Figure 11. Examples of stars identified in tlie BDF spectroscopy 
of an J?-drop M star sample, as described in the text. The absorp- 
tion features of TiO 5, CaH 2 and Na I as defined in table |3] are 
indicated by labels and / or grey shading on an inset magnification 
of the low wavelength features. 



of this line by defining a region to contain the line, and two 
adjacent regions as a proxy for the continuum level. This 

definition is shown in table [51 

As discussed by ISchiavon et al.l (|l997h , this feature is 
a sensitive indicator of surface gravity for stars cooler than 
3900K, and hence discriminates between dwarf stars with 



Band 

TiO 5 
CaH 2 

Na I EW 



SI 

7042 - 7046A 
7042 - 7046A 

8169-8171A 
8209-8211A 



W 

7126 - 7135A 
6814 - 6846A 

8172-8209A 



Table 3. The re gions used in spectroscopic indices defined by 
iReid et al. 1 1|1995^ a nd the Na I doublet equivalent width defined 
bv lFaber fc Frenc 3 l|l980l '). Indices are defined as i? = Fw/Fsi 
where the spectral region W measures the line strength and the 
sidebar 51 is indicative of the continuum strength. 



a surface gravity of \og{g) « 5 and giant stars of the same 
spectral class with a typical surface gravity log((?) — 1.0. 

Giant stars of class M4 and later are rare, and not ex- 
pected to substantially affect the number counts of faint 
stars, but the presence of a significant subpopulation of gi- 
ants could theoretically impact our distance distributions 
(which are based on distance modulus, assuming that the 
stars are dwarfs). As expected, the Na I equivalent width of 
all eighteen spectra are entirely consistent with dwarf stars 
(having EW> 5), and inconsistent with those of M giants 
(which would have an EW< 1.5). 

Given their optical colours (discussed in section |4)| the 
M dwarfs described here are expected to be comparable to 
or slightly lower in metallicity than M dwarfs observed in a 
local sample. To test this hypothesis, we consider the ratio of 
two molecular line indices. Low metallicity dwarfs have lower 
ratio s of CaH to TiO than higher metallicity dwarfs (Gizij 
allowing their ratio to be used as a crude metallicity 
indicator. 

The strength of these lines is usually gauged by a serie s 
of narrowband indices first measured bv lReid et al.l l|l995l ). 
These indices measure the strength of an absorption feature 
relative to a pseudo-continuum flux measured nearby. The 
definition of the TiO 5 and CaH 2 indices used here is given 
in table |3] 

On figure [12] we plot the ratio of line indices measured 
on the 18 spectra of the combined GOODS and BDF sam- 
ples. There is little evidence for significantly sub-solar metal- 
licities (i.e. [Fe/H]»i-2). Around 40% of the sample (7 of 
18) are consistent with slightly sub-solar metallicities (i.e. 
[Fe/H]~-1), while the same fraction are completely consis- 
tent with the indices for a local (solar metallicity) sample. 
20% of the sample (4 stars) show anomalously strong CaH 
features given the strength of their TiO absorption bands. 
Of these, two also show Ti05< 0.3 suggesting that they are 
of spectral type M6 or later (consistent with their photo- 
metric colours). 

We caution the reader that it is possible that these spec- 
tral indices may be affected by noise due to the subtraction 
residuals of night sky emission, and by night sky absorption, 
which isn't removed by our data reduction. Spectroscopy 
aimed at identifying high redshift galaxies is routinely ob- 
tained in such a manner as to optimise night sky line sub- 
traction. By offsetting the telescope between exposures, such 
that the target moves along the slit, the frames can be com- 
bined during reduction in such a way as to produce a simul- 
taneous measurement of the sky fiux that can be subtracted 
from the object spectrum. The measured signal to noise on 
the spectra account for the error in the sky line subtrac- 
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Figure 12. The spectroscopic abundance indices of faint M- 
stars identified during spectroscopic follow-up of high redshift 
dropout samp l es. Bl ue circles indicate the GOODS sample of 
IVanzella et~aLl l|2006l ) while red triangles show the BDF sample 
presented in this work. Typical errors on these indices are better 
than 10%. The spectroscopic indexes of local M and L stars are 
shown for reference. Points indicate the local (solar meta l licity ) 
sample of iHawlev. Gizis. fc ReidI 11996') and 'Rci d et al.l lll995l) 
while small symbols indicate the sample of Gizis ( 199^ (diamonds 
— solar metallicity, crosses - [Fe/H]=-1, asterisks - [Fe/H]=-2). 



tion, and hence given our requirement of a signal to noise 
exceeding ten, we estimate statistical error on our spectral 
indices of less than 15%. However, we note that the indices 
are measured in a spectral region affected by a weak sky- 
line complex, and that our spectra are of comparatively low 
resolution. 

While we do not attempt further quantitative analysis 
of our spectra, we note that there are other spectral features 
that may be accessible even in these faint sources. In partic- 
ularly, it may be possible in future to measure the metallicity 
using the KI doublet feature at 7700 A. This doublet is ob- 
served to be extremely strong in the low meta l licity extreme 
subdwarfs of 'Burgasse r. Cruz, fc Kirkpatrickl (|2007l ). While 
the line is not as strong in early M dwarfs as in later sub- 
types, we note that our spectra appear qualitatively similar 
to those from the nearby near-solar metallicity with no sign 
of the enhancement in line strength that might be expected 
at low metallicities. 

The typical metallicity of the sample as derived from 
18 spectra is consistent with that determined from their op- 
tical colours (figure \Sjp). Given that none of the eighteen 
stars in our spectroscopic sample are of very low metallicity 
([Fe/H]<-2.), the fraction of significantly sub-solar sources 
in our photometric sample is likely to be small (< 6%), con- 
sistent with the suggestion of two sources out of 124 (< 2%) 
based on their optical colours. 



8 A METAL-RICH POPULATION OF 

M-CLASS STARS AT HALO DISTANCES 

Given the evidence of their optical and infrared colours, com- 
bined with photometry of a number of examples, it appears 
that the M star population identified at faint magnitudes 
(and hence in extragalactic surveys) is relatively metal-rich 
([Fe/H]>-1) compared to the metal-poor halo population 
that might be expected to dominate star counts at large dis- 
tance. In fact, their metallicities are more typical of those 
observed locally in the thick disk. This suggests that there 
is only a weak metallicity gradient in the class M and later 
dwarfs identified as photometric dropout sources, although 
this analysis does not exclude the possibility that a low 
metallicity M-star population also exists at faint magni- 
tudes. 

Given this relatively metal-rich population, it is possi- 
ble to explore the co mparison with brighte r M stars further. 
As demonstrated by iHawlev et al.l (|2002l ) , for dwarf stars 
with metallicities [Fe/H]>-1, i' — z' colour increases linearly 
with subclass between spectral classes M3 and M9, fiattens 
through the early L dwarfs and then increases again. As a 
result, M stars identified as u-drops (which are typically too 
blue in i' — z' to be L stars) can be crudely classified by their 
i' — z colours. We chose to use i! — z colour for this classifi- 
cation rather than — i' for two reasons: in order to minimise 
errors on the colour, and in order to provide a second colour, 
thus beginning to constrain the SED for these sources. In 
the deep optical imaging employed here, the typical pho- 
tometric errors on this colour are of order 0.05 mag, less 
than the measured variation of colours within each subclass 
in a brighter sample from the SDSS (determined by com- 
paring spectral typing based on template fitting, spectral 
indices and ph otometric colour) which is typically 0.2 ma g 
in these fihers. (|Bochanski et aLll2007l : IHawlev et all 120021 '). 
This scatter in i! — z' colour within subclasses arises from 
variation in physical properties of the stars in question. A 
continuum in the metallicity and temperature of individual 
stars will affect their line indices and photometry, leading 
to a gradual range of characteristics in a field population 
rather than the discrete set of subtypes that might be ex- 
pected from a single age, single metallicity population such 
as a globular cluster. This uncertainty in intrinsic colour 
limits photometric classification to ±1 subclass. The surface 
density of faint stars decreases rapidly with increasing M 
star subclass, since a smaller volume is probed to the same 
magnitude limit. As a result stars in the range M4-M5 are 
expected to outnumber later classes. 

We make an approximate classification of our sample 
of ii-drop M stars by assigning each to th e closest subclass 
in typ ical colour, based on the templates of lBochanski et al] 
(|2007l ). Asfi gures[T3]and[T4]iIlustrate, early M stars (classes 
M3-M4) are found spanning the entire range from the satu- 
ration limit to the faintest magnitudes at which the GOODS 
data can distinguish unresolved sources from compact galax- 
ies. 

Later M stars (spectral class M5 and later) are in- 
trinsically fainter and so our survey probes progressively 
smaller volumes with increasing subclass. These are not seen 
at bright magnitudes but become increasingly common at 
i > 22. The latest stellar class seen in our v-drop samples 
is M7. Stars of classes later than M7 are intrinsically faint 
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and so will only be detected in a relatively small volume 
probed locally. In addition, the very red i' — z colours of 
these sources is likely to lead to them dropping below the 
limit of the i! band used for selection. Later classes may be 
detected rather as i'-band dropouts. 

We note that M3 stars should formally be too blue to 
be selected given our v — i' colour cut. As discussed above, 
classification by i' — z colour is at best an approximate 
method, and the large 'M3' sample may be explained as 
a combination of M4 stars with bluer than average i' — z' 
colour and genuine M3 stars with a redder than average 
V — i' colour. Spectral classification would be necessary to 
classify individual sources. Interestingly these 'MS' sources 
are more numerous than the 'M4' sample suggesting that 
scatter in the i' — z colour is skewed to the blue rather 
than symmetrical which may be indicative of low metallicity 
as discussed in section |4] If so, a slight correction must be 
applied to the absolute magnitude of each subclass resulting 
in slightly lower estimated distances. However, while this 
will apply a systematic shift to the distribution of distance 
moduli, it will not change the distribution. 

As expected, the two probable low metallicity candidate 
ti-drops in the GOODS-N field defy classification based on 
i! — z' colour. 

Using the typic al absolute magnitude for each class 
(as determined by iBochanski et al.l |2007| ) we can esti- 
mate the distance modulus, and hence distance to each 
star, essentially using the relation between colour and abso- 
lute magnitude d etermined from the loca l disk population 
l|Bochanski et al.| [2007: Ifa wlev et al.ll2002l ') to derive a pho- 
tometric parallax for each star in the sample. While this 
has traditionally been carried out using the 1/ vs V — J 
relati on, we chose instea d to use i' vs il — z (in common 
with iHawlev et all I2OO2I ) since these sources are brighter, 
and hence their magnitudes are more reliable, in the red- 
wards bands. 

The distribution of source distances from the sun are 
shown in the right-hand panels of figures [13] and 1141 We 
are observing class M3-M4 stars at distances between 1 and 
14kpc. The distant tail of this population is well above both 
the scale height of 275 pc for the thin disk, and the 1.5 kpc 
scale height of the thick disk that ma y be more appropriate 
for disc M dwarfs (|Zheng et al.ll200lf ). 

The most extensive survey of faint M dwarfs to date 
has been a survey of high galactic latitude fiel ds with 
TJST/WFPC l|Bahcall et al.ll 19941 : IZheng et al.ll200ll ). These 
authors probed about 1400 stars, reaching a depth of /voga = 
23.7 (Jab ~ 23.2). Their finding of a sharp decline in the 
number of stars with a relat ively short scale h eight of < 
1 kpc above the galactic disc IZheng et al.l l|200ll ) seems in- 
consistent with our fiat number counts out to a heliocentric 
distance of 10 kpc. Clearly our sample is smaller by a con- 
siderable margin, but is also significantly deeper and hence 
highly complete (>98%) at large heliocentric distances. The 
fact that a similar distribution of number counts is observed 
in two widely separated fields argues against Galactic sub- 
structure as an explanation. 

A fiat number density of stars with increasing magni- 
tude is not easily reproduced given standard models of the 
galactic density profile. In order to explore this, we con- 
structed Monte Carlo simulations of a sample of 50 stars 
randomly drawn from a population of 1 x 10® stars dis- 



tributed according to a known density profile. Profiles con- 
sidere d include the thick disc density law of IZheng et al.l 
(200 ll) which takes the form of a sech^(z/h) law with a 
scale height h=300pc above the Galactic disk and 'model 
B' of iReid et al.l (|l996l ) which considers a composite model 
of thin disk and thick disc (both with sech^ density profiles 
and scale heights of 350pc and 1.5kpc respectively) together 
with a halo density profile (scaling as approximately r~"^'^ 
with radius from the galactic centre). Variations on these 
models were also considered with varied scaling between disc 
and halo components, and a variety of halo power laws. In 
each case, number counts are predicted to increase sharply 
at faint magnitudes, since the density law varies less rapidly 
than the volume element at a given magnitude. The absence 
of such an upturn suggests that this population is falling off 
more rapidly than those expected from standard disk pro- 
files, despite extending many scale heights above the plane 
of the disk, or alternately that the conventional components 
of galactic structure need to be rethought. Our results might 
be interpreted as suggest that either the scale height of the 
thick disk has been significantly underestimated (as also sug- 
gested by Rcid ct al. (1996)) or that a third distant but rela- 
tively metal-rich component is needed to explain the popula- 
tion. However, we do caution that our stellar number counts 
in the most distant bins are small and that analysis of larger 
deep surveys are needed. 

We note with interest that the presence of a large pop- 
ulation of massive astrophysical compact halo objects (MA- 
CHOs) that are faint in the optical and contribute up to 20% 
of the ha lo mass has been inferred from microlensing anal- 
yses (e.g. iTisserand et al.ll2007l : iBennettllioOsI : lAlcock et al.l 
'2000', and references therein). These sources have been in- 
ferred to have a mass of ~O.5M0, similar to the mass of 
an early to mi d M-dwarf star. A ccording to the widely 
used 'S'-model l|Alcock et al.1 12OO0I '). the MACHO objects 
occupy a spherical, isothermal halo of core radius 5 kpc. In 
at least one case, a lensing star has been confirmed photo- 
metrically and spectroscopically to be an M4-5 dwarf star, 
with a red colour similar to thos e discussed in this analysis 
(Vf555w - Ifsuw = 3.2 fVegal. lAlcock et al.llioOll '). While 
the M dwarf population studied here may explain only a 
fraction of the proposed MACHO mass contribution to the 
dark halo, we probe only a small subset of the halo M star 
population with the colour and fiux limits applied in our 
study. We note that these cool, red stars would not easily 
be identified in the relatively shallow and blue monitoring 
imagin g employed by the MA CHO (|Alcock et al.ll2o6ol ') and 
EROS (|Tisserand et al.ll2007l ') collaborations. 



Our survey also highlights both the utility and the dif- 
ficulty of relying on the relatively small deep fields designed 
for extragalactic observations. While such surveys are ex- 
tremely sensitive and benefit from multiwavelength obser- 
vations - often including coordinated observations in the in- 
frared and deep spectroscopic follow-up - they also suffer 
from being necessarily small in size. The newer generation 
of ground-based, wide-field imagers have enabled larger deep 
fields to be observed, but under seeing conditions that make 
the effective separation of stars and high redshift galaxies 
(which become comparable in number density at faint mag- 
nitudes) impossible. The variation of 40% in the M star num- 
ber counts between the GOODS fields suggest that larger 
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Figure 13. (L eft) The i' — z' colours of i>-drop sele cted low mass stars in the GOODS-N field. M-star subclass increases linearly with 
i' — z' colour l|Hawlev et al.l 12001 IWest et al.ll2005h enabling a crude classification of stars based on their photometry. (Right) The 
distance to those stars, classified into subclass es by i' — z ' colou r and applying the spectroscopic parallax calibrated typical absolute 
magnitudes for each subclass as determined by IWest et al.l (|2005l ). The maximum distance probed by our photometric sample for each 
subclass is indicated by an arrow. 
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Figure 14. (Left) The i' — z' colours of i)-drop selected low mass stars in the GOODS-S field. (Right) The distance to those stars, 
classified into subclasses by i' — z' colour and applying the typical absolute magnitudes for each subclass as determined by I West et al.l 
l|2005l ). As in figure [131 



deep fields should be observed at HST resolutions in order 
to average out structure in the inner halo or outer disc. 



Clearly further work in larger fields with both deep and 
high resolution imaging (e.g. COSMOS. fScoville et al]|2006l ') 
is required to strengthen the observational result that the 
number density of sources remains constant to faint mag- 
nitudes, which may imply a rapid cut-off in the thick disk 
density profile at large distances above the Galactic plane. 



9 CONCLUSIONS 

We have investigated the properties of a faint population of 
M stars selected using a photometric drop method similar 
to that used to identify high redshift galaxies. The M-dwarfs 
in this study are up to five magnitudes fainter than those 
considered by previous authors, and hence lie at distances 
up to ten times further than those investigated before. 

This population extends well above the galactic plane, 
with little evidence for a decline in number at faint mag- 
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nitudes/large distances. This is difficult to reconcile with 
models of known disk structure. 

The M-stars selected by this photometric method show 
similar colours to those within a kiloparsec of the sun, both 
in the optical and the infrared. Together with the evidence 
from line indices, this suggests that there are significant 
numbers of M-stars with moderate to high metallicities ex- 
tending well above the plane of the galaxy. 

Our sample of stars with v — i' > 1.7 should be com- 
plete for mid-M dwarfs of near-solar metallicity but is in- 
complete for low-metallicity M dwarfs many of which have 
bluer colours or are more distant. Determining the metal- 
licity of individual sources is difficult without exceptionally 
deep spectroscopy. If one or more of our sources are at low 
metallicity, they may be even cooler - and hence somewhat 
closer - than we believe, although our spectroscopic results 
indicate that fewer than 6% of M stars are likely to fall into 
this category. 
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